Transport through graphene double dots 
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We present Coulomb blockade measurements in a graphene double dot system. The coupling of 
the dots to the leads and between the dots can be tuned by graphene in-plane gates. The coupling 
is a non-monotonic function of the gate voltage. Using a purely capacitive model, we extract all 
relevant energy scales of the double dot system. 



The control of individual electrons and spins jl|, 0, 
0] has been achieved in semiconductor quantum dots. 
Graphene as a material systems lends itself for small and 
well controlled quantum systems 0, [Hj] with the addi- 
tional possible benefit of increased spin coherence times 
@. Controlled coupling of quantum dots is a prerequi- 
site for the envisioned implementation of spin qubits Q 
in such systems. 

Here we demonstrate Coulomb blockade of a graphene 
double dot system. With a model of purely capacitively 
coupled dots [8[ we extract all the characteristic energy 
scales of the double dot system. The tunnel barriers be- 
tween the dots and the source or drain contact, as well as 
the tunnel barrier separating the two dots can be tuned 
by graphene in-plane gates. We show that the tunnel 
coupling is a non-monotonic function of the plunger gate 
voltage as expected from transport experiments on gated 
graphene nanoribbons 0, [TO, EL 03 • 

The graphene flakes are produced by mechanical cleav- 





FIG. 1: (a) Scanning force micrograph (SFM) of the dou- 
ble dot structure studied in this work. The outline of the 
graphene regions is highlighted by the dashed lines. The dots, 
labelled by L and R, have a diameter of 90 nm and the con- 
striction between them has a lithographic width of 30 nm. 
The dots are connected by 20 nm wide constrictions to source 
and drain contacts. The gates, labelled by CL, GL, GC, GR 
and CR, are located 40 nm from the structure, (b) Model for 
the analysis of the double dot system. C m is the mutual ca- 
pacitive coupling between the dots. Cgl,l and Cgr,r are the 
capacitances coupling the electrochemical potential in the left 
and right dot to the voltages Vgl and Vgr applied to their re- 
spective gates labelled by GL and GR in Fig. (a). Cgr,l and 
Cgl,r are the capacitances between the right gate and the 
left dot and vice- versa. The capacitances Cs and Ct> describe 
the electrostatic coupling of the source and drain contacts to 
the double dot. 



ing of natural graphite flakes and deposition on a highly 
doped silicon substrate covered by 285 nm of silicon 
dioxide [l3[. Thin flakes are found by optical mi- 
croscopy. Raman spectroscopy is used to select single- 
layer flakes [H], EH . The contacts are defined by electron 
beam lithography, followed by the evaporation of Cr/Au 
(2 nm/40 nm). Finally the graphene flake is patterned 
by denning the structure with electron beam lithogra- 
phy using a 45 nm thick PMMA resist layer, followed by 
reactive ion etching based on Ar and O2 (2:1) [l6j]. 

Fig. [T](a) shows a scanning force micrograph of the 
double dot structure studied in this work. The struc- 
ture consists of two central graphene islands forming the 
dots labeled L and R in Fig. [1] (a). They are mutually 
connected by a 30 nm wide constriction. Each dot has 
a diameter of about 90 nm. The dots are connected by 
20 nm wide constrictions to source and drain contacts. 
In addition to the doped substrate which acts as a global 
back gate, there are five graphene in-plane gates allowing 
to fine-tune the structure. The gates GL and GR can be 
used to change the number of carriers in the dots, while 
the gates CL, CR and GC are used to tune the transmis- 
sion of the constrictions and the coupling between the 
dots. For all gates, no leakage currents can be detected 
for applied voltages up to ±10 V. 

All the measurements presented in this paper are 
recorded at a temperature of 1.4 K, with V applied 
between the back gate and the graphene double quan- 
tum dot circuit. 

Fig. [2] (a) shows a measurement of the current through 
the double dot as a function of the voltages applied to 
gates GR and GL for an applied bias voltage Vbias = 
500 /xV. The honeycomb pattern characteristic for the 
charge stability diagram of a double dot [8[ can be ob- 
served. Elastic transport through the double dot is only 
possible in the case where the electrochemical potentials 
in both dots are aligned mutually and with the Fermi 
energy in the leads. This is the case at the so-called 
triple points in the corners of the hexagons of constant 
charge configuration. A plot of the same measurement 
displaying the current on a logarithmic scale [Fig. [2] (b)] 
makes the connecting lines between the triple points visi- 
ble. Along these lines, only one of the dot levels is aligned 
with the Fermi energy in the leads, leading to current by 
co-tunneling processes. The current through the edges 
of the hexagons can be suppressed by changing the volt- 
ages applied to gates CL and CR in such a way that the 
barriers are less transparent. 
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FIG. 2: Current through the double dot as a function of the 
voltages Vgr and Vgl applied to the gates GR and GL for 
Vbg = V, Vol = -4.65 V, Vgc = -1.2 V and Vcr = V. 
(a),(b): Representation in linear (a) and logarithmic (b) scale 
for a bias voltage Vbias = 500 fiY applied between source and 
drain contacts, (c) Representation in linear scale for a bias 
voltage Vbias = —4 mV. (d) Illustration of the measurement 
displayed in (c), with annotation of the quantities used to 
deduce the energy scales of the system. 



When applying a higher bias voltage Fbias = — 4 mV 
between source and drain contacts, the triple points 
evolve into triangular-shaped regions of increased current 
@] [Fig. [2] (c)]. We use the model of purely capacitively 
coupled dots, presented in Fig. [1] (b), to estimate the 
energy scales of the system, assuming that the applied 
source-drain voltage drops entirely over the double dot 
system. Fig. [2](d) shows a schematic of the measurement 
from Fig. [2] (c), with indication of the different param- 
eters used in the following to estimate the energy scales 
of the system. The extension of the triangular-shaped 
regions allow the determination of the conversion factors 
between gate voltage and energy. The lever arm between 
the left gate GL and the left dot is c^gl,l = Vbi as / SVql — 
0.13 and between the right gate GR and the right dot 
<^gr,r = Vbias /5Vgr = 0.12. The lever arms between 
the left gate and the right dot and vice- versa are deter- 
mined from the slope of the co-tunneling lines delimiting 
the hexagons: c^gl,r = 0.06 and q^gr,l = 0.05. The di- 
mensions of the honeycomb cells AVgl and AVgr give 
the capacitances between the gate GL and the left dot 
Cgl — e/ AVql — 2.4 aF, and between GR and the right 
dot Cgr = e/AV^GR = 2.0 aF. The corresponding total 
capacitances of the dots are Cl = Cgl/<^l = 18.0 aF and 
Cr = Cgr/&r = 17.5 aF. This corresponds to single- 
dot charging energies Eq = <^gl,l * AVgl = 8.9 meV 
and Eq = «gr,r • AVgr = 9.2 meV, which is compara- 
ble to the values found for graphene single dots of sim- 




0.5 1 

V GC (V) 



1.5 



FIG. 3: (a),(b): Examples for two extreme cases of strong 
(a) and weak (b) mutual coupling between the dots. The 
voltages applied to the gates are Vbg = V, Vcl = —4.65 V 
and Vcr = V , V G c = -1.9 V (a) and V G c = V (b). 
(c) Coupling energy as a function of the voltage Vbc applied 
to the gate GC. The voltages applied to the other gates are 
Vbg = V, V CL = -4.65 V and V C r = 1 V. The ranges of Vgl 
and Vgr are adjusted for each value of Vbc using appropriate 
values based on the measured lever arms in order to always 
stay always at the same triple point. 



ilar size [17]. The coupling energy between both dots 
can be extracted from the splitting of the triple points: 
Eg = a GL ,L • AVql = a GR ,R ■ AVgr = 2 -5 meV. For this 
temperature and gate voltage range the tunnel coupling 
energy between the two dots is below the experimental 
resolution. 

By changing the voltage applied to the central plunger 
gate labelled GC in Fig. [T] (a), we are able to change 
the coupling between both dots. Fig. [3] (a,b) show an 
example of the charge stability diagram for two extreme 
cases of strong (a) and weak (b) coupling between the 
dots. The only difference between both measurements is 
the voltage applied to the gate GC [(a): Vgc = —1.9 V, 
(b): Vgc = V]. The coupling energy between the dots 
changes by more than a factor of two: (a): E™ = 4.2 meV 
(b): Eq =1.7 meV. The single-dot charging energies do 
not change significantly compared to the case shown in 
Fig. [21 assuming the lever arms are still the same. 

Fig. [3] (c) presents a closer analysis of how the cou- 
pling energy changes with Vgc- It displays the coupling 
energy as a function of the voltage Vgc fc> r the same triple 
point, followed through the Vgr - Vgl parameter space 
as Vgc is changed. The lever arms are assumed to be the 
same as in Fig. O The strength of the inter-dot coupling 
shows a non-monotonic behavior as a function of applied 
Vbc- the coupling energy starts at quite small values, 
increases by a factor of two for more positive Vgc , before 
decreasing again. This is in agreement with recent exper- 
iments on graphene nanoribbons 0, [TO, EH, EH , showing 
a strongly non-monotonic dependence of the current on 
gate voltage, with many sharp resonances. The conduc- 
tance has been shown to vary over orders of magnitude 
as a function of in-plane gate voltage @, 03 • Here, we 
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rather probe the electrostatic landscape between the two 
dots and the corresponding change in coupling is only 
a factor of two for a given resonance. It remains to be 
seen how the tunneling transmission of the constriction 
can be linked to the electrostatic coupling it provides be- 
tween two quantum systems. 

In conclusion, we have demonstrated Coulomb block- 
ade in a graphene double dot system. The coupling be- 
tween both dots, as well as the transmission of the con- 
strictions connecting the dots to the leads, can be tuned 
by graphene in-plane gates. We have shown that the cou- 



pling between the dots is a non-monotonic function of the 
applied gate voltage. Finally, a model of purely capaci- 
tively coupled dots allowed to extract the relevant energy 
scales of the system. The presented results may be seen 
as a promising development towards the realization of 
spin qubits in graphene. 
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